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ABSTRACT 



Within the cold dark matter (CDM) framework tidal dwarf galaxies (TDGs) cannot contain dark matter, so the recent 
results by Bournaud et al. (2007) that 3 rotating TDGs do show significant evidence for being dark matter dominated is 
inconsistent with the current concordance cosmological theory unless yet another dark matter component is postulated. 
We confirm that the TDG rotation curves are consistent with Newtonian dynamics only if either an additional dark 
matter component is postulated, or if all 3 TDGs happen to be viewed nearly edge-on, which is unlikely given the 
geometry of the tidal debris. We also find that the observed rotation curves are very naturally explained without any 
free parameters within the modified Newtonian dynamics (MOND) framework if inclinations are adopted as derived by 
Bournaud et al. We explore difi'erent inclination angles and two different assumptions about the external field effect. 
The results do not change significantly, and we conclude therefore that Newtonian dynamics has severe problems while 
MOND does exceedingly well in explaining the observed rotation curves of the 3 TDGs studied by Bournaud et al. 
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1. Introduction 

Dwarf irregular galaxies have the tendency to form dur- 
ing galaxy interactions within the extended tidal tails, and 
are therefore dubbed Tidal Dwarf Galaxies (TDG, Mirabel, 
Dottori & Lutz 1992; Due et al 2000, Braine et al. 2000). 
Tidal tails transport the angular momentum and energy 
away from the merging sub-systems which build up the 
larger galaxies we see today, and gravitational instabili- 
ties within the tidal tails lead to local collapse and star 
formation. Observations of local interacting galaxies have 
shown the formation of TDGs to be quite common, in some 
cases dozens of condensations are seen in the tidal tails 
(Weilbacher et al. 2000), with a few ones having a mass 
typical of dwarf galaxies (Braine et al 2000). Based on a 
sample of 6 local interacting galaxies, Delgado-Donate et al. 
(2003) estimate that at most a few long-lived TDGs forms 
per merger. Monreal-Ibero et al (2007) studied in detail the 
stability of external star forming regions in ULIRGs, and 
concluded that they are good TDG candidates. The rela- 
tion in terms of evolution beween TDGs and the other kinds 
of dwarf galaxies is presently a matter of debate (see, e.g.. 
Due & Mirabel 1998, Okazaki et Taniguchi 2000, Bournaud 
& Due 2006, Metz & Kroupa 2007). In particular, Okazaki 
& Taniguchi (2000) claim that conservative assumptions 
about TDG production within a hierarchical CDM struc- 
ture formation framework imply that all dE galaxies may be 
TDGs leaving little or no room for traditional dark-matter 
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filled dwarf galaxies (but see Bournaud & Due 2006). On 
the other hand, Skillman & Bender (1995) cast some doubt 
on the idea that all dwarf ellipticals originate from dwarf 
irregulars after their gas has been blown away. 

In any case, in the Cold Dark Matter (CDM) scenario, 
TDGs should be nearly devoid of dark matter (Barnes & 
Hernquist 1992), contrary to other kinds of dwarf galaxies. 

In a recent paper, Bournaud et al. (2007) (hereafter 
B07) analyse the rotation curves of 3 TDGs belonging to 
the NGC 5291 system. They find evidence for a mass dis- 
crepancy that is unexpected within the CDM framework, 
and they put forward the hypothesis of baryonic dark mat- 
ter to explain the observations. Here we investigate the 
possibility that the observations can be explained within 
the framework of MOND (Modified Newtonian Dynamics, 
Milgrom 1983), without the need for baryonic dark matter. 

In general, disk galaxies' rotation curves do not de- 
crease in the outer parts as would be expected from the 
visible matter distribution. In the standard picture, this is 
explained by dark matter halos around galaxies. However, 
MOND is an alternative explanation where our understand- 
ing of gravity (or inertia) is changed, rather than our un- 
derstanding of the matter content of the Universe. Milgrom 
(1983) postulated that for gravitational accelerations below 
oo « 10~^cms~^ the effective gravitational attraction ap- 
proaches {gNioY^'^ where is the usual Newtonian grav- 
itational field. MOND leads to remarkable fits of galactic 
kinematics over 5 decades in mass (Sanders & McGaugh 
2002) ranging from tiny dwarfs (e.g.. Gentile et al. 2007) 
through our own Milky Way (Famaey & Binney 2005; 
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Famaey, Bruneton & Zhao 2007) and early-type spirals 
(e.g., Sanders & Noordermeer 2007), to massive ellipticals 
(Milgrom & Sanders 2003), without resorting to galactic 
dark matter. However, possible problems for MOND on 
galactic scales include the kinematics of Local Group dwarf 
spheroidals (Lokas, Mamon & Prada 2006) and galaxy 
merging timescales (Nipoti et al. 2007). 

TDGs are thus objects on galactic scales on which to 
test the MOND paradignfl Here we take into account un- 
certainties on the external field effect (EFE) and the galax- 
ies' inclinations to show that the MOND paradigm is very 
likely to correctly explain the kinematics of the TDGs. 

2. Data 

We use the HI rotation curves recently published by B07. 
They performed HI observations of the NGC 5291 sys- 
tem, using the VLA (Very Large Array) in its BnA and 
CnB configuration, yielding a spatial resolution of ^ 2.2 
X 1.6 kpc. They discuss the ring structure composed of 
collisional debris that surrounds NGC 5291, and they de- 
tect eight clumps where velocity gradients can be observed. 
Only three of these, however, are resolved enough to allow 
a more detailed investigation of their kinematics. 

The rotation curves of these three TDGs (NGC5291N, 
NGC5291S, and NGC5291SW) were derived by B07 using 
an envelope-tracing method (Sancisi & Allen 1979; Sofue 
1996; Gentile et al. 2004), which gives reliable velocity esti- 
mates for poorly resolved and/or highly inclined disk galax- 
ies. For the purpose of our analysis, we folded the two sides 
of the rotation curve Vmt (?") and the baryonic contribution 
Vha.r{r), to obtain azimuthally averaged Vrotir) and Vharir) 
curves. B07 use an inclination i of 45° for the three TDGs, 
based on the fact that their model reproduces successfully 
the morphology of the system with i — 45° , and that the 
rotation axis of the clumps differs from that of the ring by 
less than 15°. The total baryonic masses, as derived by B07, 
are: 9lJ:? x 10^ Mq for NGC 5291N, 9.3lo:9 x 10^ Mq for 
NGC 5291S, and 5 ± 1.5 x 10** Mq for NGC 5291SW. 

3. Fitting procedure 

To make a rigorous fit of these galaxies in MOND, two 
important issues to take into account are their inclination 
and the EFE of MOND (e.g., Famaey et al. 2007; Wu et al. 
2007, Angus & McGaugh 2007). The rotation curves were 
thus fitted within the MOND framework with five different 
hypotheses: 

• First, we simply applied the MOND prescription to 
the rotation curve. This approach has no explicit free pa- 
rameter: we used i, the distance and the stellar mass-to- 
light (M/L) ratio as derived by B07. The M/L ratio is not 
a crucial assumption since the stellar mass, derived from 
stellar population synthesis models, is several times lower 
than the total baryonic mass. We used oq = 1.2 x 10"^ cm 
s~^ (Begeman, Broeils & Sanders 1991). 

• Second, we left i as a free parameter. Indeed i is an im- 
portant parameter influencing the derived dynamical mass, 
and the inclination determined by B07 is only an estimate. 
The inclination we used is however constant with radius: 

^ A preliminary fit just posted on the arXiv by Milgrom (2007) 
has shown that MOND has the ability to explain the dynamics 
without resorting to disk dark matter. 



i.e., we ignored any possible warps, which are taken into 
account in the error-bars. 

• Third, we added a first estimate for the EFE, while i 
was left as a free parameter. In modified gravity (Bekenstein 
& Milgrom 1984), if one considers, as a first approximation, 
that a galaxy free-falls with a "uniform" acceleration in an 
external linear potential, the internal potential becomes po- 
larized Keplerian (Milgrom 1986, Zhao & Tian 2006, Zhao 
& Famaey 2006). Here, to simplify the treatment of the 
EFE, we follow the less rigorous approach of Famaey et 
al. (2007) in which the polarization in the direction of the 
external field (EF) as well as a term proportional to the gra- 
dient of fi are neglected. We define the internal acceleration 
as aint — Kibs/'^j where V^obs is the observed rotation veloc- 
ity and r is the galactocentric radius in the TDGs. Then, 
following the MOND hypothesis: Oint = g-[<!/ii{x), where 
5n = Vhar/^ *he Newtonian acceleration, ^{x) = x/{l+x) 
(Famaey & Binney 2005) and x — {a-mt +fl'ext)/ao (instead 
of the usual x — flint /ao), where ^cxt is the EF. We roughly 
estimate ^cxt due to NGC 5291 as follows: from the total HI 
profile AV20 (Malphrus et al. 1997), and assuming that the 
rotation curve of NGC 5291 stays flat out to the distances 
of the TDGs, we find ffext ^ (0.5AF2o/sm(i))2 /Rq, where 
Rq is the projected distance of the TDGs from the centre of 
NGC 5291. This gives t/oxt ~ O.Sap. We then assume f/oxt to 
have a typical value of 0.2flo. Hence, the galaxies are mostly 
dominated by aint until the last data point, so the various 
approaches to take the EF into account are slightly degen- 
erate, which we illustrate with x = max(aint/aoj ffext/flo) 
as a second parametrization. Our two choices of x for the 
EFE bracket the range of possible effects of the EF on the 
rotation curve, the first one giving a larger EFE than the 
second. A more rigourous treatment of the effect of EF on 
rotation curves has already been done in Wu et al. (2007). 
It was found that the internal iso-potential contours are 
twisted due to the general orientations of the EF, but there 
are typically less than 10% difference among different ori- 
entations (hence the modulus of vector sums) of the EF 
and internal field. Rotation curves would become normal 
baryon-only Newtonian RCs if the system accelerates faster 
and faster in a stronger and stronger EF. 

• Fourth, we fitted the rotation curves with Newtonian 
gravity and no dark matter, as would be expected in CDM, 
still keeping i as a free parameter. 

• Finally, we fixed i at the most favourable inclination 
for the previous fits, and we derived the baryonic distribu- 
tion required to fit the rotation curves in Newtonian gravity. 

4. Results 

In Fig. [1] we show the results of the fits. Overall MOND 
gives good fits to the rotation curves of these systems. In 
particular. Fig. [1] (panels a - c) shows that with no ex- 
plicit free parameter, MOND can explain well the observed 
kinematics. As could be expected a priori, leaving j as a 
free parameter improves the quality of the fits (see Fig. [H 
panels d - f). The best-fit values are not far from the first 
assumption of i = 45°: we find i = 50^7^ degrees for NGC 

5291N, i = 34±4 degrees for NGC 5291S, and i = 39t" de- 
grees for NGC 5291SW: their average and scatter is similar 
to B07. It is interesting to note that all three TDGs fall on 
the baryonic TuUy-Fisher (BTF) relation (McGaugh 2005): 
the agreement is perfect for i as a free parameter and only 
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Fig. 1. Rotation curve data (full circles) of the 3 tidal dwarf galaxies (Bournaud et al. 2007). The lower (red) curves are 
the Newtonian contribution Vbar of the baryons (and its uncertainty, indicated as dotted lines) . The upper (black) curves 
are the MOND prediction and its uncertainty (dotted lines). The top panels have as an implicit assumption (following 
Bournaud et al.) an inclination angle of 45 degrees. In the middle panels the inclination is a free parameter, and the 
bottom panels show the fits made with the first estimate for the external field effect (EFE) (Section [S]). 



slightly worse for i = 45° (see Fig. [2|). TDGs add to a long 
list of classes of galactic objects that satisfy the BTF rela- 
tion. It is puzzling that an empirical relation works so well 
for galaxies of a wide range of formation histories unless it 
is the direct prediction of a realistic law of gravity relating 
to the instantaneous baryon distribution, as MOND is. 

The fits with the two EFE estimates (see Fig. [TJ panels 
g - i, for the first estimate) give good fits, whose best-fit 
parameters are displayed in Table[TJ The first EFE estimate 
gives slightly worse fits than the second; we thus conclude 
that the EFE due to NGC 5291 in the 3 TDGs cannot be 
too large: it has to be smaller than ^ (0.1 — 0.2)ao. 

We realise that the EFE is a function of distance Rq 
from the host galaxy NGC 5291. The 3 TDGs are at pro- 
jected distances ranging from 58 to 75 kpc; given the small 
range of projected distances, and since the 3D distance is 
unknown, the EFE estimates are made using their average 
projected distance, 65 kpc, instead of their individual pro- 
jected distances, which results in errors on the EFE ^ 15%, 
much smaller than the overall uncertainties on the actual 
EFE value. We ran models of the EFE with variations of 
~ 15%, and found very similar results. Similarly, the differ- 
ence between the EE effects on the two sides of the TDG 
rotation curves, is smaller than the uncertainty on its value, 



so we ignored it and worked on the azimuthally averaged 
rotation curve. 

In Fig. [1] the curves with EFE show a Keplerian fall-off 
just outside the present data, a signal which is falsifiable 
if future observations could extend the kinematics of these 
TDGs. This is not due to tidal truncation, but to the EFE. 
A rotating disk inclined with an angle 9ef with respect 
to the EE will generally precess around the axis of the EF 
(Wu et al. 2007, Famaey et al. 2007). The precession means 
that these TDG disks in MOND are not expected to be ex- 
actly aligned with the same inclinations, justifying our fits 
where i is left as a free parameter. We note that continued 
accretion from the tidal matter and torquing of the TDGs 
may, however, affect the shape of the rotation curve at large 
radii, so the solid lines at large radii only show the case for 
unperturbed satellites. 

We note that large uncertainties still exist in the mass 
distributions and inclinations of these TDGs, and thus in 
their dynamics. The scenario proposed by B07 within the 
CDM paradigm leaves some liberty as to the exact geome- 
try and encounter. In particular, the inclination of the TDG 
on the sky planes could be more edge-on, and the inferred 
enclosed mass required to explain the observed velocities 
smaller, decreasing the need for dark matter or modified 
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Table 1. Best-fit values of the 5 fits for each TDG. i is the 
inclination, and the EFE estimates refer to Section [31 The 
1-cr uncertainties in the fitted incUnations are based on the 
statistics. 
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Fig. 2. Baryonic TuUy-Fisher relation (baryonic mass vs. 
circular velocity). The small full circles are the disk galaxy 
data from McGaugh (2005). The 3 TDGs studied here are 
shown as blue empty circles (i = 45°) and red stars (i as a 
free parameter). 



gravity. The overall quality of Newtonian, purely baryonic 
fits with i — 90° and with the baryonic distribution taken 
from B07 (not shown in Fig. [1]) is worse than the MOND 
fits with a free inclination, but it is only slightly worse than 
the MOND fits with i = 45°: in NGC5291N the MOND 
fit is much better, in NGC5291S it is a bit worse, and in 
NGC5291 the two fits are equivalent. However, given the 
geometry of the NGC 5291 system, it seems unlikely (but 
possible) that all 3 TDGs studied have such high inclina- 
tions. Also, the radial distribution of baryons, essentially 
the HI gas, is also poorly determined, and can be chosen 
to fit the observed rotation curves in Newtonian gravity, 
instead of being selected a priori. We have made such a 
fit (Fig.ini), by selecting exponential profiles for the surface 
density of the interstellar gas distributed in a disk. The en- 
closed gas mass within each radius is also plotted, together 
with the gas surface density, compatible with the observa- 
tions. The deduced rotation curves are consistent with the 
data, within the error bars. Let us note that the fit is still 
rough, no effort was made to fit more components (such as 
the stellar component and the molecular gas) in addition 
to the HI gas. 

Given the (relatively) poor quality of the rotation 
curves, a more rigorous study of the EFE and of the pos- 
sible warps is difficult, since they are not well constrained 
observationally. The uncertainty on the EFE is unlikely to 
significantly affect our conclusions given that the EFE and 
the no-EFE results do not differ much. 



5. Conclusions 

We analysed the rotation curves of 3 TDGs presented by 
Bournaud et al. (2007) and we found that MOND explains 
their kinematics very naturally. Within the MOND frame- 
work, no additional unseen matter is required in these 3 
TDGs. The dark matter required in Newtonian dynamics 
implies a large amount of unseen baryonic matter, which is 
very unconventional in the CDM context. 

We explored different possibilities for the inclination an- 
gles and the external field effect (EFE) due to NGC 5291. 



The results do not change significantly: obviously the fits 
are improved with one extra free parameter (the inclina- 
tion), and we find the external field to be smaller than 

(0.1 — 0.2)ao, within the a priori expectations from the 
kinematics of NGC 5291. Therefore, MOND can fit the 
kinematics of the 3 TDGs presented by Bournaud et al. 
(2007) for inclinations similar to those that they estimated 
and for realistic values of the external field. 

We also found that in the (rather unlikely) case where 
all 3 TDGs are seen edge-on, Newtonian fits with no dark 
matter are possible but slightly worse than the MOND fits 
with i = 45°, and not as good as the MOND fits with i as 
a free parameter. This Letter has thus shown that TDGs 
pose severe challenges to the current standard cosmolog- 
ical theory and that very recent observations of rotating 
TDGs are consistent with a non-Newtonian interpretation. 
TDGs thus hold important clues on fundamental physics 
and therewith require further research. 

Acknowledgements. We thank S. McGaugh and G. Angus for useful 
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